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SUMMARY 


As aircraft become larger and lighter due to design requirements for 
increased payload and improved fuel efficiency, they may also become much 
more flexible, for highly flexible vehicles, the handling qualities may not be 
accurately predicted by conventional methods. This study applies two analysis 
methods to a family of flexible aircraft in order to investigate how and when 
structural (especially dynamic aeroelastic) effects affect the dynamic 
characteristics of aircraft. The first type of analysis is an open-loop modal 
analysis technique. This method considers the effect of modal residue 
magnitudes on determining vehicle handling qualities. The second method is a 
pilot-in-the-loop analysis procedure that considers several closed-loop system 
characteristics. Both analyses indicated how dynamic aeroelastic effects can 
cause a degradation in vehicle tracking performance, based on the evaluation of 
•wine simulation results. 

This report is divided into two volumes. Volume I consists f the 
development and application of the two analysis methods described above. 
Volume II consists of the presentrtion of the state variable models of the 
flexible aircraft configurations used in the analysis applications, mode shape 
plots for the structural modes, numerical results from the modal analysis, 
frequency response plots from the pilot-in-the-loop analysis and a list'.ig of the 
modal analysis computer program. 
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Appendix A. 3 
Vehicle Configurations 
and 

Modal Analysis Results 


The vehicle configurations are presented in state variable form using the 


following state and output definitions. 


x v T = [<*,%, u, 0 R , rh, t/ 2 ,rj 3 , r/ 4 , fj„ i 2 , t} 3 , rj 4 ] , 

(A.3.1' 

y T = [u, 0 R , 0 R , Tf, 0 T , Of, n z ] . 

(A.3.2) 


- where all angles are expressed in radians, all angular rates in radians per 
second, ihe forward velocity perturbation, u, has units of feet per second and 
.the generalized elastic deflections, rj , are dimensionless. The matrices which 
describe the vehicle configurations are defined in the following equations. 

i, = [AVIs. + (Bis, + |D|S, 


(A.3.3) 


X = |C]S, + |E|s,, + |F|s, 

The vectors, x p and x g , are defined Chapter 4. 

The results of the modal analysis consist of several parameters which are 
presented in matrix form. The equations which follow present the definitions 
of these modal analysis results. 

The similarity transformation which was used to scale the eigenvectors, as 
described in Appendix A.l, has the following definition. 
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T - SCALING MATRIX 

= diag[l, 1 , 1 , ft v ft 2 , 

u o 

ft 3 > ft 41 ^ |> 0* 2 > ^ 3 ) ^ 4 ] (A. 3. 4) 

The sealed system is therefore described by the following set of equations. 

x v = Tx v (A.3.5) 


i v = T(AVJT 1 i v + T[B) x p + T[D)x< 

= [AV-SCALED] i v 

+ [B-SCALEDJ jl p + [D-SCALEDJ & g (A.3.6) 

jl = [C)T >i v + [E]x.p + (FJ*, 

= [C-SCALEDJ i v + [E]2 l p + [FJ x* (A.3.7) 

The pilot lag states and gust states are defined in the following equations. 

ip = |AP] ip + [Gp]r/p (A.3.8) 

i g = (AG]i g + [G g ]»j g (A.3.9) 

The complete aircraft system state model used in the modal analysis is defined 
in the following equations. 

2L T = | i p T ' X g T ’ 2L v T 1 (A. 3. 10) 
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x = 


(AP) 0 0 

0 (AG) 0 

[B-SCALED] [D-SCALED] [AV-SCALED] 



G P 


0 

+ 

0 

% + 

G s 


0 


0 


X. 


= [AA] x + [GP] % + (GG] n g (A.3.11) 

JL= [(E) [F] (C-SCALED]]i 

= [CC]* (A.3.12) 

The modal controllability, disturbability and observability matrices are defined 


in the following equations. 

|EE] = modal matrix of [AA] (A.3.13) 

The modal controllability matrix, 

[CON] ^ [EE] 1 (GP| . (A.3.14) 

The modal disturbability matrix, 

(DIST] = [EE] 1 [GG] . (A.3.15) 

And, the modal observability matrix, 

[OBS] = [CC] [EE] . (A.3.16) 


The residues for the complete aircraft model are defined as follows. The 
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residues of the vehicle due to the ith control input are described by, 

(C-IMP RES]i ^ |[OBS] [CONJi] 7 (A.3.17) 

- where [CON] ; is the diagonal matrix formed by the ith column of the 
modal controllability matrix and [OBS] is the modal observability matrix. 

Each kj element of [C-IMP RES]; corresponds to the residue of the kth mode 
for the jth output due to an impulse in the ith control input. Similarly, for the 
disturbance residues, 

[D-IMP RESJi ^ |[OBS] [DISTJiJ 7 (A.3.18) 


Here, [DISTlj is the diagonal matrix formed by the ith column of the modal 
disturbability matrix and so each kj element of [D-IMP RESj; corresponds to 
the residue of the kth mode for jth output due to an impulse in the ith 
disturbance input. 

The remainder of Appendix A.3 contains the numerical results of the 
modal analysis and the vehicle matrices for eight data base configurations as 
described above. 
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Appendix A.4 

Aeroelastic-Structural Mode Shapes 


The mode shapes of the data base configurations are defined to be 
deflections in feet normalized so that for each mode the deflection at the 
nose (i.e. F.S. 0) is one foot. The deflections are defined relative to a sot of 
bod> axes defined in the following way. The x-axis lies along the longitudinal 
axis of the vehicle and is directed forward. The y-axis lies in the plane of the 
wings, is orthogonal to the x-a\is and points out the right wing. The z-axis is 
orthogonal to the xy-plane and is directed downward. The origin of the xyz- 
coordinate system is defined to be the center of gravity of the aircraft. 

The deflection definitions and their sign conventions that are used in this 
study are defined in f igure A. 4.1. The mode shapes for the four longitudinal 
modes used in this study are presented in Figures A. 4. 2 through A. 4. 5. 
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MOD JO SHAPE SKIN CONVENTIONS 



Mode Shape 

/ 



<t> x ^ mode deflection (feet) 


<t >' , mode slope (feet/feet) 
i 1 X ^ generalized deflection (dimensionless) 


\ 


Figure A. 4. 1 

Mode Shape Sign Conventions 




•oos are normalised to 


* one foot displacement 


»t the nose 


, , F %ureA.4.2 

AeroeJastic-Structural Mode 


s &ape * Mode 1 
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deflection (ft) 
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MODE 2 



Figure A.4.3 

Aeroeltttic-Structur&I Mode Shape * Mode 2 





deflection (ft) 
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MODE 4 



* deflections are normalized to a one foot displacement at the ncse 


Figure A.4.5 

Aeroelastic-Structural Mode Shape - Mode 4 


251 


Appendix \.5 

OCM Frequency Responses 
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Figure A.5.1 (eon’t) 

OCM Frequency Responses - Config. 1 
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Figure A. 5. 2 (con’t) 

OCM Frequency Responses - Config. 2 
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Figure A.5.2 (con’t) 

OCM F requency Responses • Cocfig. 2 
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Figure A.5.3 (coo t) 

OCM Frequency Responses - Config. 3 
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Figure A. 5.3 (con’t) 

OCM Frequency Responses • Confi^ 3 
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Figure A> '< (con’t) 

OCM Frequency Responses - Config. 3 
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Figure A.5.4 

OCM Frequency Responses - Config. 4 
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Figure A.5.4 (eon’t) 

OCM Frequency Responses - Config. 4 
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Figure A.5.4 (con’t) 

OCM F requency Responses - ConBg. 4 
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Figure A.5 4 (con’t) 

OCM Frequency Responses - Config. 4 
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Figure A.5.5 (con’t) 

OCM F requency Responses - Config. 5 
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Figure A S 6 (eon’t) 

OCM Frequency Respc’.ses - Config. 6 
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Figure A. 5. 6 (con’t) 

OCM Frequency Responses - Config. 6 
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Figure A..V7 

OCM f requency R*»ponsej - Config. 7 
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Figure A.5.7 (con’t) 

OCM Frequency Responses - Config. 7 
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Appendix A.O 

Listing of Modal Analysis Program 


The modal analysis computer program which was used for the open-loop 
analysis method is presented here. This program is written in FORTRAN and 
makes use of a modular structure to simplify its use and to allow the reader to 
more easily follow the procedure. Documentation is integrated with the code 
and a variable name listing and description is included to further aide the 
reader. 

Kven though the program was written to be flexible, it is based on a linear 
system analysis package of subroutines that is unique to the Purdue computer 
system. The package was developed at Purdue and is entitled LSLIB. 
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LIST OF UARIABLES USED IN MODAL ANALYSIS PROGRAM 




VARIABLE NAME t VARIABLE TYPE T DESCRIPTION 

UTILITY UARIABLE U5ED IN SUBROUTINE PHAG 
AUGMENTED SYSTEM MATRIX 
GUST EQUATION SYSTEM MATRIX 


A 

AA 

AG 

AMASE 


AP 

ARRAY 

AS 

AU 

B 

B 

BB 

BS 

C 

CASE 

CC 

CNM 

CON 

CS 

D 

BIST 


T 

r 

T 

T 

T 

T 

T 

T 

T 

T 

T 

T 

t 

T 

T 

t 

T 

T 

r 

T 

T 

t 

T 

t 

t 

T 

T 

t 

T 

T 

T 

T 

T 

T 

t 

T 

T 

T 

T 

t 

T 

t 

t 

T 

T 

T 

t 

T 

T 


REAL 

REAL ARRAY 


T 

T 

T 

t 


REAL ARRAY 
REAL ARRAY 


REAL ARRAY 
REAL ARRAY 
REAL ARRAY 
REAL ARRAY 
REAL ARRAY 
REAL 

REAL ARRAY 


REAL ARRAY 
REAL ARRAY 


REAL ARRAY 
REAL ARRAY 


REAL ARRAY 


COMPLEX ARRAY 


REAL ARRAY 
REAL ARRAY 
COMPLEX ARRAY 


ARRAY USED IN FORMING MAG/PHASE FORM OF 
COMPLEX MATRICES 

PILOT EQUATION SYSTEM MATRIX 

UTILITY ARRAY USED IN SUBROUTINE PHAC 

UTILITY ARRAY USEO IN SCALING MATRIX AU 

VEHICLE EQUATION SYS'SEM MATRIX 

VEHICLE EQUATION CONTROL MATRIX 

UTILITY UARIABLE USED IN SUBROUTINE PHAC 

AUGMENTED CONTROL MATRIX - USED ONLY IF 
THE PILOT CHARACTERISTICS ARE NOT BEING 
MODELED IN THE ANALYSIS 

UTILITY ARRAY USED IN SCALING MATRIX 1 

OUTPUT MATRIX ASSOCIATED WITH THE 
VEHICLE STATES 

ALPHA -NUMERIC USED TO TITLE OUTPUT 

AUGMENTED OUTPUT MATRIX FORMED BY 
COMBINING C. fc\ AND F 

UTILITY VECTOR USED TO DETERMINE ELEMENTS 
OF C ASSOCIATED WITH PLUNGE ACCEl ORATION CN2) 

CONTROLABILITY MATRIX FORMED FROM THE MODAL 
MATRiX (EE) AND THE SYSTEM CONTRQ; 

MATRIX (BB OR CP) 

UTILITY ARRAY USED IN SCALING MAiRIX C 

VEHICLE EQUATION DISTURBANCE MATT < 

DI3TURBA3ILITY MATRIX FORMF" F :OV THE MODAL 
MATRIX (EE) AND THE S STEM WOE 

MATRIX (GG) 
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05 

T REAL ARRAY 

t 


t 

T 

E 

T REAL ARRAY 

T 


T 

T 

EE 

T COMPLEX ARRAY 

T 


T 

T 


t 

T 

El 

T COMPLEX ARRAY 

T 


T 

? 

EIC 

t COMPLEX ARRAY 

T 


T 

T 

EIGMAT 

T COMPLEX ARRAY 

T 


t 

T 

ElhU 

T COMPLEX ARRAY 

T 


T 

T 

EMM 

T REAL ARRAY 

T 


T 

T 


T 

t 

F 

t REAL ARRAY 

T 


T 

T 

FIELD 

t REAL 

T 


T 

T 


t 

T 

FIELDC 

T REAL 

T 


t 

T 


T 

T 

FMM 

T REAL ARRAY 

T 


t 

T 


T 

T 

FQRMAA 

t SUBROUTINE 

T 


T 

t 

FORMBB 

T SUBROUTINE 

T 


t 

T 

FORfICC 

T SUBROUTINE 

T 


T 

T 

FORMCA 

T SUBROUTINE 

T 


T 

T 


T 

T 

FORMNZ 

T SUBROUTINE 

T 


t 

t 


t 

T 

CG 

T REAL ARRAY 

T 


T 

T 

GGA 

T REAL ARRAY 

T 


T 

T 


T 

t 

GP 

r REAL ARRAY 

T 


T 

T 

CPA 

T REAL ARRAY 

t 


1 

T 

GRAU 

t REAL 

T 


T 

t 

H BLANK 

f REAL 

t 


t 

T 


r 

t 

I 

T INTECER 

T 


T 

? 

IMPMAl 

t COMPLEX ARRAY 

T 


t 

r 


T 

t 


UTILITY ARRAY USED IN SCALING MATRIX D 

OUTPUT MATRIX ASSOCIATED UITH THE CONTROLS 

MATRIX OF EIGENUECTORS GF MATRIX AA (ALSO 
CALLED THE MODAL MATRIX) 

IMUERSE OF MATRIX EE (INUERSE MODAL MATRIX) 

UECTOR OF EIGENUALUES OF MATRIX AA 

DIAGONAL MATRIX OF EIGENUALUES F.IC 

INUERSE OF MATRIX EIGMAT 

UTILITY UECTOR USED TO DETERMINE ELEMENTS 
OF E ASSOCIATED UITH PLUNGE ACCELERATION (N2) 

OUTPUT MATRIX ASSOCIATED UITH THE DISTURBANCES 

USED IN PRINTING REAL MATRICES (ONLY USED IN 
MPRINT. A LSLIB ROL n INE) 

USED IN PRINTING COMPLEX MATRICES (ONLY USED 
IN CMPRNT > A LSLIB ROUTINE) 

UTILITY UECTOR USED TO DETERMINE ELEMENTS 
OF F ASSOCIATED UITH PLUNGE ACCELERATION (NZ) 

FORMS AUGMENTED SYSTEM MATRIX AA 

FORMS AUGMENTED CONTROL MATRIX BB UHEN NP«0 

FORMS AUGMENTED OUTPUT MATRIX CC 

FORMS AUGMENTED CONTROL MATRIX GPA ANO 
AUGMENTED DISTURBANCE MATRIX GGA 

FORMS PLUNGE ACCELERATION ELEMENTS OF THE 
OUTPUT MATRICES C. E AND F 

DISTURBANCE EQUATION UHITE NOISE MATRIX 

AUGMENTED SYSTEM UHITE NOISE DISTURBANCE 
MATRIX 

PILOT EQUATION UHITE NOISE MATRIX 

AUGMENTED SYSTEM UHITE NOISE CONTROL MATRIX 

GRAUI TATIONAL ACCELERATION 

USED IN PRINTING MATRICES (ONLY USED IN 
MPRINT AMD CMPRNT. LSLIB ROUTINES) 

DO LOOP COUNTING UARIABLE 

UTILITY ARRAY FOR OUTPUTTING RESIDUES 
(IMPMAT - IMPULSE RESIDUE MATRIX) 
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IP 

IT TILE 

J 

JP 

JTITLE 

K 

KP 

L 

LIME 

LX 

M 

MATOUT 

MODAL 

N 


INTEGER 

IMTEGER 

IMTEGER 

IMTEGER 

IMTEGER 

IMTEGER 

IMTEGER 

IMTEGER 

IMTEGER 


T 
T 
T 
T 
t 
T 
T 
T 
T 
T 
t 

r 

T 
t 
T 
T 
T 
T 
T 

T REAL 
t 

T IMTEGER 

t 

t 

T COMPLEX hRRAY 

T 

T 

t SUBROUT I ME 

t 

t 

t IMTEGER 
t 


DO LOOP COUNTING UftR I ABLE 
ALPHA-NUMERIC USED TO TITLE MATRICES 
DO LOOP CGUMY IMG UARIABLE 
DO LOOP COUMTIMG UARIABLE 
ALPHA-MUMERIC USED TO TITLE MATRICES 
DO LOOP COUMTIMG UARIABLE 
DO LOOP COUMTIMG UARIABLE 
DO LOOP COUMTIMG UARIABLE 


T 

T UARIABLE USED IM HEADER* A LSLIB ROUTIME 
T USED TO TITLE qUTPUT 
T 

T DISTANCE FROM NOSE TO C.C. OF AIRCRAFT (FT) 

T 

t # OF COLUMNS IM MATRIX ARRAY IM 
T SUBROUTINE PHAC 
T 

T UTILITY MATRIX FOR OUTPUT TING COMPLEX 
t MATRICES III MHG/PHASE FORM 

T 

T DETERMINES MODAL MATRIX. INUERSE MODAL 
t MATRIX AMP £ IGEMUALUES OF A MATRIX 
T 

T H OF COLUMNS IN MATRIX MATOUT IM 
T SUBROUTINE PHAG 



t 


t 





MACC 

T 

IMTEGER 

t # OF AG COLUMNS 




T 


T 





HAGR 

T 

IMTECER 

t # OF AC ROUS 





T 


T 





NAME AC 

T 

IMTEGER 

T 

ALPHA-MUMERIC 

NAMIMC 

THE 

AC MATRIX 


t 


t 





MAMEAP 

t 

IMTEGER 

T 

ALPHA-NUMERIC 

MAMIMC 

THE 

AP MATRIX 


T 


t 





NAMEAU 

T 

IMTEGER 

T 

ALPHA-MUMERIC 

MAM I MG 

THE 

AU MATRIX 


T 


t 





NAMEB 

t 

IMTEGER 

t 

ALPHT -NUMERIC 

MAMIMC 

THE 

B MATRIX 


t 


T 





NANEC 

t 

INTEGER 

t 

AtPHA-NUMERIC 

MAMIMC 

THE 

C MATRIX 


t 


t 





MAMED 

t 

IMTEGER 

t 

ALPHA-MUMERIC 

MAMIMC 

THE 

D MATRIX 


T 


t 





MAMEE 

T 

IMTECER 

T 

ALPHA-MUMERIC 

MAMIMC 

THE 

E MATRIX 


t 


t 





NAMEF 

t 

IMTECER 

T 

ALPHA-MUMERIC 

MAMIMC 

THE 

F MATRIX 


T 


T 





MAMECC 

t 

IMTECER 

T 

ALPHA-NUMERIC 

MAMIMC 

THE 

CC MATRIX 


T 


T 





MAMECP 

t 

IMTECER 

T 

ALPHA-MUMERIC 

NfyilMG 

THE 

CP MATRIX 


T 


t 





MAPC 

t 

IMTECER 

T 

* OF AP COLUMNS 




t 


t 

H OF AP ROUS 




MAPR 

t 

IMTECER 

t 
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t 

t 

KWC 

t INTEGER 

t 4? OF AU COLUMNS 


t 

T 

MAUR 

T INTEGER 

T 4 OF AU ROWS 


t 

T 

mbc 

t INTEGER 

T # OF B COLUMNS 


t 

T 

NBR 

T INTEGER 

T # OF B ROWS 


t 

T 

MC 

T INTEGER 

T £ OF CONTROLS IN UEHICLE EQUATION 


T 

T 

MCC 

T INTEGER 

T OF C COLUMNS 


T 

T 

MCDU 

T INTEGER 

T NC ♦ ND «■ NU 


T 

T 

MCP 

T INTEGER 

T USED IN PRINTING MATRICES (ONLY USED IN 


T 

T MPRINT AND CMPRNT » LSLIB ROUTINES) 


t 

T 

HCPC 

t INTEGER 

T USED IN PRINTING MATRICES (ONLY USED IN 


T 

T MPRINT AND CMPRNT# LSLIB ROUTINES) 


T 

T 

NCR 

t INTEGER 

t #GF C ROWS 


T 

T 

ND 

T INTEGER 

T It OF DISTURBANCES IN UEHICLE EQUATION 


T 

T 

NDC 

T INTEGER 

T |t OF D COLUMNS 


T 

T 

NDR 

T INTEGER 

T It OF D ROWS 


T 

T 

NDU 

T INTEGER 

t ND ♦ NU 


T 

t 

NF' 

T INTEGER 

T It OF E COLUMNS 


T 

t 

her 

T INTEGER 

t if OF E ROWS 


t 

T 

NFC 

T INTEGER 

T # OF F COLUMNS 


T 

T 

NFR 

T INTEGER 

t If OF F ROWS 


1 

t 

NGGC 

T INTEGER 

T ft OF CG COLUMNS 


t 

t 

NCGR 

t INTEGER 

T It OF GG ROWS 


f 

T 

NGPC 

t Integer 

T It OF GP COLUMNS 


T 

T 

NGPR 

T INTEGER 

T 4 OF GP ROUS 


t 

t 

NM 

t INTEGER 

t ft OF OUTPUTS IN OUTPUT EQUATION 


t 

T 

f<P 

t INTEGER 

T FLAG FOR USING PILOT EQUATION IN 


t 

t THE ANALYSIS <NP«OJ NO PILOT - NP*NCi PILOT) 


t 

T 

NPDU 

t INTEGER 

T NP ♦ NO ♦ NU 


t 

T 

NRA 

t INTEGER 

t ROW DIMENSION OF MATRIX MATOUT IN 


T 

T SUBROUTINE PHAG 


t 

T 

NRB 

T INTEGER 

T ft OF RIGID BODY STATES IN UEHICLE EQUATION 


T 

t 

NROU 

T INTEGER 

T ROW DIMENSION OF MATRIX ARRAY IN 
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r t 

t t 

MS/ t INTECER T 

r t 

t t 

MU T INTEGER T 

t t 

MW T INTEGER T 

r r 

T T 

t T 

r;wx T INTEGER r 

t r 

MWXY T IMTECER T 

* t 

MWY T IMTECER T 

T T 

MX t IMTECER T 

T T 

T T 

T f 

MY t IMTECER T 

T T 

t f 

t T 

MZ T IMTECER T 

T T 

t T 

T T 

T T 

OBS T COMPLEX ARRAY T 

T T 

T T 

PHAG T SUBROUTINE T 

T T 

f f 

PHI T REAL ARRAY t 

1 T 

PH IP t REAL ARRAY T 

T T 

RESID T SUBROUTINE T 

T T 

T T 

t T 

SCALE T SUBROUTINE 

t r 

T T 

T T 

STEP T SUBROUT t 
T T 

T T 

STP r COMPLY T 

T T 

T T 

3TPMAT t COMPLEX AkkmV t 
T r 

T T 

SUBT T REAL ARRAY T 

T f 

T T 


SUBROUTINE PHAC 

ft OF STRUCTURAL MODE STATES MODFLED IN 
THE UEHICLE EQUATION 

#0 F STATES IM THE UEHICLE EQUATION 

ROW DIMENSION OF UEHICLE SYSTEM MATRICES 
(THIS IS THE LARGEST POSSIBLE MODEL UMlCH 
AM BE ANALYZED WI THOUT ALTERING THE PROGRAM) 

NW ♦ NX 

NW ♦ MX ♦ NY 

NW ♦ NY 

ROW DIMENSION Or' PILOT SY51EM MATRICES 
(THIS IS THE LARGEST PILOT MODEL WHICH 

the program can handle without alteration) 

ROW DIMENSION OF DISTURBANCE SYSTEM MATRICES 
(THIS IS THE LARGEST GUST MODEL WHICH THE 

PROGRAM CAN HANDLE WITHOUT ALTERATION) 

• 

ROW DIMEN 3 ION CF THE OUTPUT MATRICES 
(THIS IS THE LARCEST POSSIBLE NUMBER OF 
OUTPUTS WHICH THE PROGRAM CAM HANDLE 
WITHOUT ALTERATirN) 

Q8SERUA6ILITV MP T RIX FORMED FROM THE MODAL 
MATRIX (EE) AND irtC OUTPUT MATRIX (CO 

OUTPUTS A COMPLEX MATRIX IN PHASE AND 
MHCNITUDE FORM (PHASE IN DECREES) 

MODE SHAPES FOR EACH OF THE STRUCTURAL MODES 

MODES SLOPES FOR EAC“ OF THE STRUCTURAL MODES 

DETERMINES THE IMPULSE RESIDUE MP T RIX OF A 
SYSTEM GIUEM ME CONTROLABIL I TY OR 
DISTURBA9ILITY AND CBSERUMBILITY MATRICES 

SCALES THE MATRICES OF A SYSTEM (IE. DYNAMIC. 
CONTROL. DISTURBANCE AND OUTPUT MATRICES) 

WITH A DIACUNAL SCALING Matrix 

FORMS STEP RESIDUE MATRIX CIUEM IMPULSE 
RESIDUE MATRIX AND L'ECTOS* UF EICEMUALUES 

UTILITY ARRAY FOR FORMING RESIDUE MATRICES 
(USED IN SUBROUTINE STEP) 

UTILITY ARRAY FOR OUTPUTTING RESIDUES 
(STPMAT » STEP RESIDUE MATRIX) 

UARIABLE USED IN PRINTING THE HEADING OH 
THE OUTPU T 


i 
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T t REAL ARRAY 

t 

TERM1 t REAL 
t 

TERM2 t REAL 
t 

TINU t REAL ARRAY 

t 

TITLE T REAL ARRAY 
t 
T 

UZERQ T REAL 

T 

WORK T SUBROUT I hE 

T 
T 
T 

WK T REAL ARRAY 

t 

X T REAL ARRAY 

t 
T 


t SCALING MATRIX FOR SCALING THE UEHICLE STATES 
T 

T INTERMEDIATE VARIABLE IN SUBROUTINE PHAG 
t 

T INTERMEDIATE VARIABLE IN SUBROUTINE PHAG 
T 

T INVERSE OF MATRIX T 
T 

T VARIABLE USED FOR PRINTING THE HEADING ON 
T THE OUTPUT 
T 

T CRUISE VELOCITY OF AIRCRAFT 
T 

T DETERMINES THE CONTRGLADILITY* DI5TUR3A3ILITY 
T AND OBSERVABILITY MATRICES OF A SYSTEM GIVEN 
T THE SYSTEM MATRICES AND THE MODAL MATRIX 
T 

T UTILITY VECTOR USED IN SEVERAL SUBROUTINES 

r 

T TRANSFORMATION VECTOR USED IN DETERMINING 
T ELEMENTS OF THE C MATRIX ASSOCIATED WITH 
T THE PLUNGE ACCELERATION (NZ) 
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t 

ANALYSIS OF FLEXIBLE AIRCRAFT DYNAMICS USING MODAL ANALYSIS METHODS T 


======= T 


T 

THIS PROGRAM USES A MODAL ANALYSIS TECHNIQUE TO DETERMINE THE T 

IMPULSE RESIDUES OF A STATE UAPIABLE MATH MODEL OF A FLEXIBLE t 

AIRCRAFT. IT DOES THIS ANALYSIS BY FORMING AN AUGMENTED SYSTEM T 

FROM THE MATRIX EQUATIONS DESCRIBING THE UEH1CLE DYNAMICS AND T 

THE FRECUENCY CHARACTERISTICS OF A HUMAN PILOT AND ATMOSPHERIC T 

DISTURBANCES. T 

T 

1 

T 

UEHICLE MATRIX EQUATION - T 

t 

X-DGT = [ AUJ*X + tBl»U + CDI*H T 

T 

y = cci*X + CE]*U ♦ EF3* U t 

T 

T 

PILOT MATRIX EGUATION - T 

T 

U-DOT = C AP]»U + C GPI^ETA t 

T 

T 

DISTURBANCE MATRIX EQUATION - f 

T 

U-DOT = CAG1*U + CGG1*ETA T 

r 

T 


/////////.' 


PROGRAM ANAFLEX( INPUT, OUTPUT. TAPE5=INPUT. TAPE6=GUTPUT. 
• TAPE1 , TAPE2, TAPES) 


: : : : : DIMENSION STATEMENTS : : s : s 


DIMENSION 

DIMENSION 

AU(20»20}» 

CQ5.20), 

BC20.5) . 
E( 15.5). 

0(20.5) 

F(I5»5) 

DIMENSION 

AP(5.5). 

CP(5,5) 


DIMENSION 

AG(5.5), 

GG(5,5) 


DIMENSION 

DIMENSION 

AA(30» 30). 
GPA(30>5). 

BB(25,5)t 
GGA(30» 5) 

CC( 15.30) 

COMPLEX 

COMPLEX 

EIG<30), 

C0N<30.5>, 

EE(30,30). 
0BS( IS, 30). 

EK30.30) 

DI3K30.5) 


C 

COMPLEX IMPMAT ( 30, 15). STPMAT<30. 15) 
T(20, 20). 


C 


DIMENSION 


TINUC20, 20) 
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c 


c 


c 


c 


c 


c 

c 


c 

c 


c 

c 

c 

c — 

C m • m ■ 
* « • 

c 

c 


c 

c««» 

c*#* 

c 

c 

c 


DIMENSIGN 

DIMENSION 

DIMENSION 

TITLEC7). 
PHI (4)» 

X( 1 * 20) * 

CASE(4) * 
PHIPC4) 
CNM( 1*20)* 

SUBT(2) 

ENMC1.5). FNM( 1*5) 

DIMENSION 

AS(20* 20 ) » 

BS(20* 5) * 

CS( 15* 20) * DS(20*5) 

DIMENSION 

DIMENSION 

AMASE(30* 2)* 
ICNAMEC2) * 

I TITLE (2 ) . 
IDNAMEC2) * 

JTITLEC2) 

I0NAME(2) 

COMPLEX 

COMPLEX 

MAT0UT(30* 30) 
EIGMAT(30> 30 ) 

* EINUC30* 

30) 

DIMENSION 

DIMENSION 

DIMENSION 

AR(20*2Q) * 
BR(20* 5) * 
CR( 15*20)* 

ART (20* 20 ) * 
BT( 20* 5) * 

CT ( 15* 20 ) 

ATRC20* 20) * AT (20 *20) 
DR(20* 5) * DT(2J0*5) 

LOGICAL 

TEST (20) 



DIMENSION 

DIMENSION 

A3(20* 20) • 
C3( IE 20>* 

B3(20* 5) * 
E3( 15*5)* 

D3(20*5) 
F3( 15*5) 

DIMENSION 

ATK20.20) 



DIMENSION 

DIMENSION 

(20.20)* 
Ci ;15.20)» 

BN(20*5). 
EN( 15*5)* 

DMC20.5)* 
Ff1( 15*5) 

REAL 

LX 



EXTERNAL 

PHAG 




COMMON BLOCKS ::::: 

THESE VARIABLES ARE USED IM PRINTING THE HEADINGS ON THE OUTPUT 

COHMON/MPRI NT' NCP. FIELD. NCPC . F I ELDC . HBLANK 
COMMON/HEAD/ TITLE. CASE. LINE 

THESE VARIABLES ARE THE ACTUAL WORKING MATRIX DIMENSIONS 
USED IN MOST OF THE SUBROUTINES 

COMMON/DIM/ NU. NC. ND. NM. NP 


C: : s s 
C 

c*** 

C 

c 

c. 

c*** 

c*»» 

c 

c 

c 

£ 

c 

c 

c 


SET ROW DIMENSIONS ::::: 


THESE VARIABLES ARE THE ROU DIMENSIONS OF THE SYSTEM 

ROU DIMENSION * SIZE OF MATRIX SET IN DIMENSION 

THESE ROU DIMENSIONS MUST BE AT LEAST' AS LARGE 
AS THE WORKING DIMENSIONS 


AU(NW.NW), 

AP(NX.NX). 

AGCNY.NY). 

C(NZ.NX). 

AA(NWXY.NWXY). 

GPA(NWXY.NX). 


B(NU.NX). 

GP(NX.NX) 

GC(NY.NY) 

E(NZ.NX). 

CD ( NWX. NX). 
CGA(NUXY.NY) 


D(MU.NY) 


F(NZ.NY) 

CC(NZ.NWXY) 


MATRICES 

STATEMENTS 
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C 

MM 3 20 
NX 3 5 
NY 3 S 
NZ 3 15 
NUY 3 NM ♦ NY 
NMXY * NU ♦ NX ♦ NY 
C 

c t 

C PHYSICAL PARAMETERS PARTICULAR TO THIS PROBLEM T 




C*S==S=*5S: 




C NUMBER OP STRUCTURAL MODES MODELED 

NSM 3 4 

C STRUCTURAL MODE ::: SHAPES (PT) ::: SLOPES T 

C 


PHI ( 1 ) = .349 

$ 

PHIP(l) 

= -.02640 

PHI(2) = .400 

$ 

PHIPC2) 

= -.03300 

PHIC3) 3 .180 

$ 

PHIP(3) 

= -.03456 

PHI<4) 3 .140 

$ 

PHIPC4) 

= -.03636 


C CRUISE UELOCITY (RVSEC) * t 

UZERO 3 349.0 

C result LOCATION (FT) T 

LX 3 63.4 


C CRAUITATIQNAL ACCELERATION IFT/SECySEC) T 

GRAU 3 32.17 



C 

C 

C 

Css?** READ TITLES s s s * * 

C 

C*«* TITLE ✓ CASE ✓ SUBTITLE 
C 

REAO(S.SOl) TITLE 
WRITEC3.S01) TITLE 

501 F0RMAT(7A10) 

READ(5.502) CASE 
URITE<3.502) CASE 

502 F0RMAT(4A10) 

READ (5. 503) SUBT 
WRITE! 3. 503) SUBT 

503 F0RMAT(2A10) 

C 

C»** THIS SUBROUTINE PRINTS THE HEADER ON THE OUTPUT 
C 

CALL HEADER (SUBT) 

C 

C 

Cs * : t s READ MATRICES Its:* 

C 

(>*• NAME ✓ DIMENSION (ROU.COLUMN) / ELEMENTS 
C 


READ(S*S04) NAMEAU 



URITEC3.504) NAMEAU 
READCS.*) NAUR.NAUC 
WRITEC3. *) NAUR.NAUC 
DO 1 1=1. NAUR 

READCS.*) (AU( I, J), J=l. NAUC) 

1 URITE(3* *) CAUCI,J),J=1,NAUC) 
READ (5. 504) NAME AG 

READCS. *) NAGR.NAGC 
DO 2 I=l.NAGR 

2 READC5.*) (AG(I.J). J=1.NAGC) 
READCS. 504) MAMEB 
UR1TEC3. 504) NAMEB 

READC5. *) NBR.NBC 
HRITEC3, ») NBR.NBC 
DO 3 1=1. NBR 

READCS.*) CBCI.J).J-*1.NBC) 

3 WRITEC3. *) CBCI.J), J=1.NBC) 
READCS. 504) NAMEC 
WR1TEC3.S04) NAMEC 
READC5.*) NCR.NCC 
URITEC3. *) HCR.NCC 

DO 4 1=1. NCR 

READCS,*) CCCI, J), J=1»NCC) 

4 URITEC3, *) CCCI, J). J=l.NCC) 
READCS. 504) NAMED 
URITEC3.504) NAMED 
READCS. *) NOR.NDC 
URITEC3. *) MDR.NDC 

DO 5 1=1. NDR 

READC5. *) (DCI. J). J=1,NDC> 

5 URITCC3, *) CDCI. J), J=1,NDC) 
READCS. 504) NAMEE 
WRITEC3.504) NAMEE 
READCS,*) NER.NEC 
URITEC3. *) NER.NEC 

DO 6 1=1, NER 

READC5. *) CECI, J). J=l,NEC) 

G WRITEC3,*) CECI, J),J=1, NEC) 
READCS. 504) NAMEE 
WR I TEC 3. 504) NAMEE 
READ C 5.*) NER, NFC 
URITEC3, *) NFR.NEC 
DO 7 1=1. NER 

READC5. •) CECI. J),J=1, NEC) 

7 WR1TEC3. •) CECI. J),J=1, NEC) 
READCS. 504) NAMEGG 

504 FORMAT CA4) 

READCS,*) NGGR.NGGC 
DO 8 1=1, NGGR 

8 READCS,*) CGGCI, J), J=1,NGGC) 
C 

READC5.S04) NAMEAP 
RE AD (5.*) NAPR.NAPC 
IFCNAPR.EQ.O) GO TO 10 
C 

DO 11 1=1, NAPR 

il READCS.*) (APCI, J),J=1,NAPC) 
PEADCS.504) NAMEGP 
READC5. ») NGPR.NCPC 
00 12 l«l,NCPR 
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AU<NU,NU). 
C(NM>NU)f 
AP(NC» NC)> 
AG(ND> ND) * 
AA(NCDU.NCDU), 


SCMUtMC)# 
E(NM,NC)f 
GP(NC>NC) 
GC(ND. 1) 
BBCNDU# NC) » 


DCMUi ND) 
F(NM.ND) 


CCCNH.NCDU) 


12 READC5* *0 (GP( I> J)? J=l# MGPC) 

C 

Cs::s: SET WORKING DIMENSIONS ::::: 

C 

C«*»* THESE UARIABLES ARE THE ACTUAL DIMENSION OF THE SYSTEM MATRICES 
C 
C 
C 
C 
C 
C 
C 
10 


c 

c 

c»** 

c»«* 

c 

c 


NU = 

NAUC 

NC *-= 

NBC 

ND = 

NDC 

NM * 

NCR 

NDU = 

= NU + ND 


NCDU * NU + ND ♦ NC 


INCLUDE PILOT MODEL IN ANALYSIS - NP = NC T 

IGNORE PILOT IN ANALYSIS - NP = 0 t 


?IP = NC 

NPDU = NU ♦ ND + NP 


C " 

C THIS SECTION OF THE PROGRAM IS PROBLEM DEPENDENT 


C 
C 
C- 
C 

c- 
C:s:si 
C 

c 

13 
C 


IT IS DUE TO THE CHOICE OF USING PLUNGE 
ACCELERATION AS A SYSTEM OUTPUT UARIABLE 


FORM UECTOR-X FOR HZ EQUATION ::::: 


T 

T 

'f 

T 

T 


DO 13 1*1, NU 
XC1. I) = 0.0 

X(l.l) » -UZERO/GRAU 
XC 1.2) * LX'GRAU 
X( 1,4) = UZERO/GRAU 
XC1.9) * -PHKD/GRAU 
X(l. 10) = -PHI(2)/GRAU 
X(l, 11) * -PHI(3)/GRAU 
X(l, 12) * -PHI(4)/CRAU 
C 

C::sir COMPLETE C MATRIX UITH NZ EQUATION tSISJ 

C 

C 

CALL FOPMNZI AU, B, C. D, E. F.NU.NX, NY, NZ, X. CNM. ENM.FNM) 

C 

Z///// ••////////// ////////////////////////////////////////// 

c 

c 
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C 

CSSSSS OUTPUT SYSTEM MATRICES ::::: 

c 

C 

CALL MPR I NT ( NW . AU , MU . MU , G , S , MAMEAU. 1 ) 
CALL MPRINTCNY, AG, ND» ND. ND, G, NAMEAG. 1) 
CALL MPRINTCNU.B.NU.NC.NC.G.NAMEB. 1) 
CALL MPR IN T (NZ. C.NM.NU. 6, G.NAMEC. 1 ’ 
CALL MPRINT(NU.D»NU,ND.ND.6.NAMEB. I) 
CALL MPRINT ( NZ. E. NM. NC» NC . S. MAMEE# 1 ) 
CALL MPRINT (NZ. F.NM.ND.ND. S. MAMEF » I ) 

IF (NP.EQ. 0) GO TO 14 
CALL MPR I NT ( NX » AP . MC . NC . NC . B . NAMEAP . I ) 
14 CONTINUE 
C 

C 

cs:::5 FORM SCALING MATRIX s : i t * 

C 

C 


C*— — 

C THE T MATRIX IS PROBLEM DEPENDENT 


C 

c»«» 

c*** 

c 

c 

c 

C 

o** 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 


THE ACTUAL SCALING MATRIX, T, IS PROBLEM DEPENDENT 
BUT THE SCALING PROCESS IS ALWAYS RECOMMENDED 

NUMBER OF RIGID BODY STATES 

NRB = NU - 2.0 * NSM 


IN THIS CASE. THE SCALING MATRIX, T, IS BASED ON MODE SLOPES 


1.0 


T = DIAG 


1.0 

1/U 

1.0 

Y/XU) 

Y/XC2) 

Y/XC3) 

Y/X(4) 

Y/XU) 

Y/XO) 

Y/XC3) 

Y/X<4) 


- WHERE Y/X(I) IS THE MODE SLLi'E ASSOCIATED WITH 
THE ITH MODE OF THE STRUCTURE 

-AND- 

- WHERE U IS SOME CHARACTERISTIC UELOCITY 
( IN THIS CASE U - U2ER0 ) 


DO 15 1*1, NU 
DO 16 J=1«NU 
T(I.J) » 0.0 


16 
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15 TCI* I) * 1.0 

C 

T(3»3) = 1.0 / UZERO 
C 

DO 17 1*1. NSM 
TCI+NRB.I+NRB) = PHIP(I) 

1 7 T ( I+NRB+NSM, I+NRB+NSM) = PHIP(I) 

C 

C/ ////////////////////////////// 

c 

c 

C::::: SCALE UEHICLE SYSTEM STATES ::::: 

C 

C 

CALL SCALE ( AU, B. C. D. NU. NX, NY. NZ.T. TINU. AS. BS. CS, DS) 

C 

C 

C:s*:: CUTPvJT SCALED MATRICES ::::: 

C 

c 

CALL MPRINT(NU.T,NU.NU.G.G.*SCALINC*.;) 

C 

CALL MPRINTt NU. AU. NU. NU, 6. 6. *AU- SCALED*. 1 ) 

CALL MPRINTtNU, B. NU, NC. NC. 6. *B-SCALED*. l ) 

CALL MPRIMT(NZ»C, NM.NU.S.S, *C-SCALED*» l) 

CALL MPRINTtNU. D. NU, ND. 6. 6, XD-SCALED*. 1 ) 

C 

C 

CS : s s s SET UP TEST ARRAY mil 

C 

C 

c»*» nr * THE NUMBER OF RETAINED STATES 
C**» NT * THE NUMBER OF TRUNCATED STATES 
C 

NR * 12 
NT * MU - NR 
C 

C**» INITIALIZE THE TEST ARRAY TO .TRUE. 

C 

DO 29 1=1. NU 
29 TEST(I)= .TRUE. 

C 

IF CNT.EO.O) CO TO 42 
C 

O*** IF TEST ( I ) - .FALSET. THE ITH STATE IS 

C**» TRUNCATED IN fHE RES I DUAL I ZAT ION 

C 

TESTC6)* .FALSE. 

TEST(7)* .FALSE. 

TEST(8)» .FALSE. 

TEST(IO)* .FALSE. 

TEST(ll)* .FALSE. 

TEST<12)« .FALSE. 

C 

C 

Cm 1 1 APPLY STATE MODEL RESIDUALIZATION mil 
C 

c 

CALL SORT ( AU» B. C, D. AR. ART » ATR. A f » BR. BT. DR. DT » CR* CT , TEST . 
• NU.MX.NY.NZ.NR.NT) 
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c 

Z////S/SSS OUTPUT PARTITIONED MATRICES 

CALL MPRINT ( NU. AR, NR. NR. NR. G. *ARj«, +1 ) 

CALL MPRI NT ( NU. ART. NR. NT. NT. 6. *ART*» *1 ) 

CALL MPRINTINU, ATR. NT, NR. NR. 6. *ATR*. +1 ) 

CALL MPRINT(NU. AT.NT.NT.NT»6»*AT*»+1) 

CALL MPRINT(NW, BR. NR. NC. NC, G. *BR*, +1 ) 

CALL MPRINT(NU» BT, NT. NC. NC, 6. *BT*, +1 ) 

CALL MPRINT(NU,DR.NR,ND,ND.G.*DR*.+1> 

CALL MPRINTCNU, DT, NT, ND»ND»G»*DT*. +1) 

CALL MPRINTCNZ, CR, NM, NR, NR. 6. *CR<. +1 ) 

CALL nPRINT(NZ»CT.NM,NT»MT»6»*CTX»+l) 

C 

CALL MRR(AR, ART, ATR, AT. BR. BT. DR, DT. CR. CT, E. F. AM, BM. CM, DM, EM, FM. 

* NU»NX»NY»NZ»NR»HT»A3.B3»C3»D3,£3»F3»ATI) 

C 

GO TO 41 
C 

C*** IF NO TRUNCATION SET THE UEHICLE MATRICES TO THE MODIFIED MATRICES 
C 

42 CALL MEQCAU. AM, NU, NU, NU.NU) 

CALL MEQ(B, BM.NU.NC.NU.NW) 

CALL MCQ(C.CM,NM,NU,NZ.NZ) 

CALL MEQ(D, DM. NU.ND, NU.NU J 
CALL MEQ(E.EM.NM.NC.NZ.NZ) 

CALL MEQ(F.FM.NM.ND.NZ.MZ) 

C 

41 CONTINUE 

C 

C*** RESET NU SO THAT THE UEHICLE MODEL IS 

C*** REPRESENTED BY THE SMALLER. RESIDUALIZED MODEL 

C 

NU = NR 

NDU = ND ♦ NU 

NPDU = NP * ND ♦ NU 


C 

C 

Cssss: FORM AUGMENTED WHITE NOISE MATRICES 

C C 0 1 C GP I 

C GGA = t GG J CPA * ( 0 1 

C t 0 J (01 

C 

C IF NP » 0 : 

C 

C 

C GGA = [ GG ) GPA - t 0 J 

C (01 

C 

C 


CALL FORMGAIGG.GP.NUXY.HX.NY.GGA.GPA) 

C 

c 

Csss:: OUTPUT MOISE MATRICES sssss 

C 

c 

CALL MPRIMKNWY.GGA. NPDU. 1. 1.6. MAMEGG, 1 ) 

IF (NP.EO.O) CO TO 39 

CALL MPRINT(NUXY,GPA,UPDU.NC.MC.6.NAMEGP. 1) 
CONTINUE 


39 
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C 

c 

C: : : : : FORM AUGMENTED SYSTEM MATRIX 

c 

C 

C C AP T 0 T 0 ] 

C AA-tO TAGTO ] 

C C 8 T D t AU 1 

C 

C IF MP = 0 : 

c 

C AA * C AG t 0 ] 

C C D T AU 1 

C 

c 

r 


CALL FORMAAI AM, AG# AP# BM. DM. AA# MU. MX. MY# MUXY) 
C 

Z/// / / / */ / OUTPUT MATRIX AA 

CALL MPRINTCNUXY»AA.NPDU.NPDU,8,6.*AA*, 1) 

c 

IP (NP.HE.O) CO TO 18 
C 

C 

C::::s FORM AUGMENTED 3 MATRIX SJJJJ 

C 

C ONLY IF NP - 0 S 
C 

C B3 = I 0 I 
C t B 1 

C 

C 

c 

CALL FORMBBT BM. BB* MU. MX. MUY) 

C 

Z///////// OUTPUT MATRIX BB 

CALL ilFRINTtNHX, BB. NDU, NC, 6. 6. *3B*. 1 ) 


C 

c 

C::::: FORM AUGMENTED C MATRIX 

c 

C CC * t E t F t C 1 
C 

C IF NP » 0 : 

C 

C CC =* t F T C I 

C 

C 

c 

18 CALL FORNCCICM.LM.FM.CC.NW.NX.NY.NZ.NWXY) 

C 

Z///////// OUTPUT MATRIX CC 

CALL MPR INT ( HZ . CC . Ntl . NPnU . 5 . G . XCCX . t ) 


C 



CSS*:: DETERMINE E-UALUES AMO E-UECTORS OF MATRIX AA 

C . 

C*** THIS SUBROUTINE USES GEIGEN UHILH OUTPUTS THE r-U'»H* , g 
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c 

CALL MODALCAA.EIC. EE. El. MUXY) 

C 

Z////////y OUTPUT E-UECTOR MATRIX 

CALL CMPRMT (NUXY - EE. MPDU. MPDU. 6.6. *EE*» 1 ) 


C 

C- 

r 

C 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

21 

L 

22 

c 


DETERMIME COMTROLABILITY. DISTLRBABILITY AMD 
OBSERUABTL ITY MATRICES 


CCM = MI MU * B 
DIST = MIMU * D 
OBS = C * M 

-WHERE M = MODAL MATRIX 


If (NP.EQ.O) CO TO 21 

CALL WORK (CPA. CC. GCA. EE. El. MX. MY. NZ.MWXY.COM. OBS, DIST) 
CO TO 22 

CALL WORK (BB.CC.CCA.EE. El. MX, MY.MZ.riWXY.COM. OBS. DIST) 
CONTINUE 


Zs /////////// OUTPUT COM. DIST, OBS MATRICES 
C 

DATA ICNAME/* CONTROL^* LABILITY */ 

DATA IDNAME/* DISTURB*. ^ABILITY *✓ 

DATA IOMAME/* OBSERUA*. *B ILITY */ 

CALL CMPRNT ( NUXY. COM, MPDU, MC. G, 6. *COM*. 1 ) 

CALL PHAG( CON. MUXY. MPDU, NC. ICMAME.AMASE.NATOUT.NWXY) 
CALL CMPRMT ( NUXY, DIST. MPDU. 1 . G. 6, *DIST*. 1 ) 

CALL PH ACC DIST. MUXY. MPDU. I . IDMAME. AMASE. MATOUT.NWXY) 
COLL CMPRMT ( MZ. OBS. MM. MPDU , 6. G * *OBS <. 1) 

CALL PHAGCOBS.MZ.MM.NHDU. I ONAME . AMAGC • M A TOUT , MZ ) 

C 

c 


C:::s: DETERMIME MODAL IMPULSE AMD STEP RESIDUE MATRICES 5 

c 

C 

C RESMAT - OBS > < COM (OR DIST) 

C 

C 

C 


tit 
• ♦ • 


G01 


C 

c*«» 

c 

c 


URITE(S.SOl) 

FORMAT(//v. 

* 10X,*«» RESIDUE MATRICES 

A 10X. X* *•••*»■ 

A 10X.**»» RUUS CORRESPOND . Q MODES *•**,/, 

A 10X.**»« COLUMNS CORRESPOND TO OUTPUTS 

- 10X.*** **•••*) 


I * OF UHITE MOISE CONTROL IMPUTS 
DO 19 1=1. NC 

CALL PE3ID(C0N.033« N2. NX. NUXY. IMPMAT. I) 



301 


c 

C///'/ OUTPUT IMPULSE RESIDUE MATRIX 
C 

URITECB.B02) I 

602 FORMAT C.'.SX.P *, /. 

-> 5X» P IMPULSE RESIDUE MATRIX FOR CONTROL - *. 13.', 

c " * f) 

CALL CMPRNTCMUXY,INPMAT.r?»DU.NM.S.E,PC-lMP RES P.l) 

C 

DATA ITIfLE/-* CONTROL^* X RESIDUE *' 

CALL PHAGC IMPMAT » NUXY.NPDU. NM. ITITLE.AMASE.MATOUT.MUXY) 

C 

CALL STEP ( I MPMAT . E I G. NUXY . M2 . STPMAT . E I CMAT . E I MU ) 

C 

C.'A/v/ OUTPUT STEP RESIDUE MATRIX 
C 

URITE(6.£04) I 

S04 FORMAT (✓/✓. 5X, P 

- 5X.P STEP RESIDUE MATRIX FOR COTTROL - P.I3./. 

* 5 X.p *,✓) 

C 

CALL CMPRNT (NUXY. STPMAT . NPDU, MM. E. B. PC-STP RES X. 1) 

C 

CALL PHAGC SI PMAT • NUXY. MPDU. MM, ITITLE.AMASE.MATOUT.MUXY) 

C 

ts cormrac 
c 

C»** I * * OF UHITE MOISE DISTURBAMCES 
C 

I - 1 

C 

CALL RESIDCDIST.OBS.MZ.MY.MUXY. IMPMAT. I) 

C 

C///// OUTPUT IMPULSE RESIDUE MATRIX 
C 

URITECB.B03) I 

U03 FORMAT (/'•'•', 5X, P 

» .'.5X.P IMPULSE RESIDUE MATRIX FC? DISTURBANCE - 

A SX.P P,/) 

C 

CALL CMPRNT ( MUXY. IMPMAT. MPDU. NM.G.S. PD- IMP RES P. 1 ) 

C 

DATA JTITLE-'PDI STURBAMCp , PE RESIDUE */ 

CALL PHAGC IMPMAT. NUXY. MPDU. MM. JTI TLC . AMASE . MA TOUT. MUXY ) 

C 

CALL STEPC IMPMAT ,£IC. MUXY.M2, STPMAT. El CMAT . EIMU) 

C //✓✓/ OUTPUT STEP RESIDUE MATRIX 
C 

URITEC6.S05) I 

GCC FORMATC///.SX.P *,,, 

- 5X.P STEP RESIDUE MATRIX FOR DISTURBANCE - P.T3.'. 

^ » 5X.P - / ) 

CALL CNPRNTC NUXY. STPMAT. MPDU. MM. 6. E.PD-STP PES P. 1 ) 

c 

CALL PHAGC STPMAT . MUXY. MPDU.IIM. JTITLF. AMAST . MAT 'L -iUW) 

C 

REWIND 1 
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STOP 

END 


l 
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i : : t : m 1 1 1 1 1 > < 1 1 1 : s : : t : t : : 1 1 s : : i > m t : m : 1 1 1 : s : 1 1 : 1 1 : 1 1 
SUBROUT I ME FORMAA 


THIS SUBROUTIME FORMS THE AUGMENTED MATRIX - 

I AP T 0T 01 
AA 3 I 0 T AG T Cl 
C B T D T AU] 

- FROM THE MATRICES AP. AG. AU. B AMD 0 
IF HP = HC 

AMO IE AUGMEMTED MATRIX - 

AA = I AC T 01 
I 0 t AU1 

- FROM THE MATRICES AG. AU AMD D IF MP 3 0 


CALL STATEMEMT ARGUMEMTS - 

AU ! MU X MU MATRIX OF ROU OIMEMSION MU 

AC : MO X MO MATRIX UF ROU DIMEMSIOM MY 

C AP : MC X MC MATRIX OF ROU DIMEMSIOM MX 

C B : MU X MC MATRIX OF ROU DIMEMSIOM MU 

C D : MU X MO MATRIX OF ROU OIMEMSIOM MU 

C AA : MCOU X MCDU MATRIX OF ROU DIMEMSIOM MUXY 

C 

c 

c 

SUBROUTIME FORMAA! AU. AC, AP. B, D. AA. MU. MX, MY, MUXY) 
DIMEMSIOM AU(MU.MU). AGINY.liV). AP(MX.MX) 

DIMEMSIOM B(MU.MX) ; D(MU.MY). AACMUXY, MUXY) 

C 

COMMOM /DIM/ NU. MC. MD. MM. MP 
C 

MPDU 3 MU ♦ MD * MP 
C 

C INITIALIZE MATRIX AA TO ZERO 

C 

DO 10 I 3 l , MPDU 
DO 20 J* l » MPDU 
20 AA( I. J) 3 0.0 

10 CONTINUE 

C 

IF (MP.EQ.O) CO TO 35 
C 

C FORM UPPER LEFT PARTITIOh OF MATRIX AA IF MP 3 NC 

C 

DO 30 IM.MP 
DO 40 J-l.NP 
40 AA(I.J) 3 API 1 1 J) 

30 COMTIMUE 

C 

C FORM LOUER LEFT PARTITION OF MATRIX AA IF MP ■ MC 

C 


DO SO I-ltKJ 
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no 60 J* 1 * NC 

GO AA( I+NP+ND. J) = B(I.J) 

50 CONTINUE 

c 

C FORM CENTER PARTITIION OF MATRIX AA IF NP = NC. 

C OR THE UPPER RIGHT PARTITION OF MATRIX AA IF NP = 0 

C 

35 DO 70 1=1. ND 

DO 80 J=1 » ND 

80 AA( I+NPt J+NP) = AG(I.J) 

70 CONTINUE 

C 

C FORM LONER RIGHT PARTITION OF MATRIX AA 

C 

DO 90 1=1. NU 
DO 100 J=1 * NU 

100 AA( I+NP+ND. J+NP+ND) = AUCI.J) 

90 CONTINUE 

C 

C FORM LONER CENTER PARTITION OF MATRIX AA IF NP = NC 

C GR THE LOUER LEFT PARTITION OF MATRIX AA IF NP = 0 

C 

DO 110 1=1. NU 
DO 120 J= 1 t ND 

120 AA( I^NP+ND. J+NP) = D(t.J) 

110 CONTINUE 

C 

RETURN 

END 
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SUBROUTINE* FORMBB 


THIS SUBROUTINE FORMS THE AUGMENTED MATRIX - 

BB = t 0 J 
( B 1 

- FROM THE MATRIX B AND THE ZERO PARTITION 
HAS ND ROUS 


CALL STATEMENT ARGUMENTS - 

B : NU X NC MATRIX OF ROW DIMENSION NW 
BB : NDU X NC MATRIX OK ROU DIMENSION NWY 
NX : COLUMN DIMENSION OF MATRICES B ANr BB 


SUBROUTINE FORMBBI B. BB. NU. NX. NUY) 
DIMENSION B(NU.NX). BB(NUY.NX) 

COMMON /DIM/ NU. NC. ND. NM. NP 

NDU a NU ♦ ND 

INITIALIZE MATRIX BB TO ZERO 

DO 10 1*1. NDU 
DO SO J=1.NC 
SO BB(I.J) a 0.0 
10 CONTINUE 

. FORM LOWER PARTITION OF MATRIX BB 

DO 30 1=1, NU 
DO 40 J=1.NC 

40 BB( I+ND. J) = B(I.J) 

30 CONTINUE 
C 

RETURN 

END 
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C: 

C 

C: 

C 

C 

C 

C 

C 

C 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 


c 


****** * * SUBROUT I ME *FORMCC* “ ‘ e ****** 6 ******** ’ 

• ••••••••••••••••a*************************** 

THIS SUBROUTINE FORMS THE MATRIX - 
CC = I E T F T C J 
- FROM THE MATRICES C. E AND F 


CALL STATEMENT ARGUMENTS - 


C 

E 

F 

CC 

MU 

NX 

MY 


MM X MU MATRIX OF ROU DIMENSION M2 

MM X MC MATRIX OF ROU DIMENSION M2 

MM X MD MATRIX OF ROU DIMENSION M2 

MM X MCDU MATRIX OF ROU DIMENSION NUXY 
COLUMN D I MENS I ON CF MATRIX C 
COLUMN DIMENSION OF MATRIX E 
COLUMM DIMENSION OF MATRIX F 


SUDROUT I NE FORMCC ( C . E . F . CC . MU . NX . NY . NZ , MUXY ) 

DIMENSION C(NZ.NU). E(NZ.HX). F(NZ.NY), CQ(NZ.NUXY) 


COMMON 'DIM/' NU. MC. ND. NM. NP 
C 

MPDU = NU ♦ ND + NP 
C 

C INITIALIZE MATRIX CC TO ZERO 

C 

DO 10 1=1. NM 
DO 20 J=1.NPDU 
20 CC(I.J) = 0.0 
10 CONTINUE 
C 

IF (NP.EQ.O) CO TO 25 
C 

C FORM LEFT PARTITION OF MATRIX CC IF NP * NC 

C 

DO 30 1=1. MM 
DO 40 J=l,NC 
40 CC(I.J) = E(I. J) 

30 CONTINUE 

C 

C FORM CENTER PARTITION OF MATRIX CC IF NP = NC 

C OR FORM LEFT PARTITION OF MATRIX CC IF NP = 0 

C 

25 DO 50 1=1. NN 
DO 60 J=1.ND 

GO CCd.JHIP) = F(I.J) 

50 CONTINUE 

C 

C FORM RICHT PAR TIT ION OF MATRIX CC 

C 

DC 70 1=1. NM 
DO 80 J=l.NU 

80 CC( I» J+NP+ND) = C(I.J) 
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70 CONTINUE 
C 

RETURN 

END 


v 
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SUBROUTINE FORMGA 


THIS SUBROUTINE FORMS THE MATRICES - 

t 0 I [ GP I 

GGA = C GG 1 : GPA = C 0 ] 

I 0 I CO] 

- FROM THE MATRICES CG AND GP 
IF NP = 0. WHERE THE LOWER 
ZERO PARTITIONS HAUE NU ROWS 
AND THE OTHER ZERO PARTITIONS 
ARE COMPATABLE UITH THE 
CORRESPONDING POPULATED PARTITIONS 

AND. THE MATRIX - 

GGA = t GG 1 
( 0 ] 

- IF NP = NC. WHERE THE LOWER 
ZERO PARTITION HAS NU ROWS 


CALL STATEMENT ARGUMENTS - 

GG : ND X l MATRIX OF ROW DIMENSION NV 
GP : NC X NC MATRIX OF ROW DIMENSION NX 
CCA : NPDU X 1 MATRIX OF ROW DIMENSION hWXY 
(IF NP = 0; NPDU = ND ♦ NU 
IF NP = NC? NPDU = NP + ND + NU) 

GPA : NPDU X NU MATRIX OF ROW DIMENSION NWXY 


SUBROUTINE FORMCA ( GG . GP , NWXY . NX , NY , GG A , GPA ) 

DIMENSION GCCNY, I). GP(NX.NX). GGA(NWXY.l). GPA(NUXY.NX) 

COMMON/DIM/ NU. NC. ND. NM. NP 

NPDU = NP ♦ NO *• NU 

INITIALIZE MATRICES GGA AND GPA TO ZERO 

DO 1 1=1. NPDU 
GGA(I.l) = 0.0 
DO 2 J= 1 . NC 
GPA(I.J) = 0.0 
CONTINUE 

FORM UPPER PARTITION OF MATRIX CCA IF MP = 0 
OF FORM CENTER PARTITION OF MATRIX GGA IF NP = NC 

DO 3 tOl.ND 
3 GGA (K HIP) = GGOO 

C 


IF (NP.EQ.O) RETURN 
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FORM UPPER PARTITION OF MATRIX CPA 

DO 4 L=l.NC 
DO 5 M=1»NC 
5 GPA(L.M) = CP(L.M) 

4 CONTINUE 

C 

RETURN 

END 
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C: : 


SUBROUTINE FORUMS 


::::::::::::::::::: 


THIS SUBROUTINE DETERMINES THE ROU OF ELEMENTS 
OF THE MATRICES C. E AND F UHICH CORRESPOND 
TO TERMS OF AN OUTPUT EQUATION WHICH DESCRIBES 
AN OUTPUT UARIABLE WHICH IS PROPORTIONAL TO 
THE FIRST DERIUATIUE OF THE STATES. UHERE- 

Y * CX + EU + FU. 


THE ROW WHICH CORRESPONDS TO THE THIS RESPONSE 
IS ASSUMED TO BE THE LAST ROW OF THE MATRICES. 

THE OPERATION OF DETERMINING THE X-DOT 
PROPORTIONAL ELEMENTS CONSISTS 

OF EXTRACTING CERTAIN ELEMENTS OF THE MATRICES A. 

B AND D WHICH APPEAR IN THE OUTPUT EQUATION. 

IN THE SPECIFIC PROBLEM DONE IN THIS STUDY. 

THE X-DOT PROPORTIONAL OUTPUT IS PLUNGE ACCELERATION. 
THE ACTUAL N 2 EQUATION CAN BE FOUND IN THE TEXT 
DESCRIBING THIS SUBROUTINE. 


CALL STATEMENT ARGUMENTS - 

A S NU X NU MATRIX OF ROU DIMENSION NW 

B : NU X NC MATRIX OF ROW DIMENSION NW 

C : NM X NU MATRIX O* - ROW DIMENSION NZ 

D S NU X ND MATRIX OF ROW DIMENSION NW 

E : NM X NC MATRIX OF ROW DIMENSION NZ 

F s NN X ND MATRIX OF ROW DIMENSION NZ 

NX J COLUMN DIMENSION OF MATRICES B AND E 
NY : COLUMN DIMENSION OF MATRICES D AND F 
X : 1 X NU TRANSFORMATION ROW UECTOR OF 
COLUMN DIMENSION NW 

CNM*. l X NM UTILITY ROW UECTOR OF NZ ELEMENTS 
FOR THE MATRIX C OF COLUMN DIMENSION NU 
ENMs 1 X NC UTILITY ROW UECTOR OF NZ ELEMENTS 
FOR THE MATRIX E OF COLUMN DIMENSION NX 
FNM! l X ND UTII T TY ROW UECTOR OF NZ ELEMENTS 
FOR THE M. RIX F OF COLUMN DIMENSION NY 


SUBROUTINE FORMNZ(A,B.C.D.E.F.NU.NX.NY.NZ, X.CNM.ENM. FNM) 
DIMENSION A(NW.NU) » B(NW.NX). C(NZ.NW) 

DIMENSION D(NW.NY). ECNZ.NX). F(NZ.NY) 

DIMENSION XC1.NW). CNM(l.NW). ENM(l.NX). FNM(l.NY) 

COMMON /DIM/ NU. NC. ND. NM. NP 


SUBROUTINE MULT ( A. B. AB. L. M, N. NA. NB. NAB ) 

COMPUTES THE MATRIX PRODUCT: AB - A * B 

A : L X M MATRIX WITH ROW DIMENSION NA 
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SUBROUTINE MODAL 

THIS SUBROUTINE DETERMINES THE EIGENUECTORS. 
MODAL MATRIX AND INUERSE MODAL MATRIX OF A 
GENERAL MATRIX. 


CALL STATEMENT ARGUMENTS - 

AA : NPDU X NPDU MATRIX OR ROW DIMENSION NUXY 
EIG : NPDU UECTOR OF EIGENVALUES OF MATRIX AA OF 
DIMENSION NUXY 

EE : NPDU X NPDU MATRIX OF EIGENUECTORS OF AA 
OF ROW DIMENSION NWXY 

El S NPDU X NPDU MATRIX INUERSE OF EE OF RGU 
DIMENSION NWXY 


SUBROUTINE MODAL (AA. EIG. EE. El* NUXY) 

DIMENSION AA( NWXY. NWXY). UK (1000) 

COMPLEX EIG(NWXY). EE (NUXY. NUXY ) » El (NWXY. NWXY) 

CCMMON /DIM*' NU. NC, ND. NM. NP 
NPDU = NU ♦ NP ND 


SUBROUTINE GEIGEN ( NA. AA. N. JOB. UK. KPR. FROM. 

EIG. NE. EE. El. NAME ) 

COMPUTES EIGENVALUES AND. OPTIONALLY. MODAL MATRIX OR 
NORMALIZED MODAL MATRIX AND INUERSE MODAL MATRIX OF A 
REAL GENERAL MATRIX 

INPUT ARGUMENTS- 

AA(N.N): REAL GENERAL MATRIX UITH ROU 
DIMENSION NA 
N : ORDER OF AA 
JOB : OPTION CONTROL 

= C : EIGENVALUES ONLY 
= 1 S EIGENVALUES AND MODAL MATRIX 
= 2 S EICENUALUES# MODAL MATRIX AMD 
INUERSE MODAL MATRIX 
= 3 S SAME AS 1 UITH NORMALIZED 
MODAL MATRIX 

= 4 : SAME AS 2 UITH NORMALIZED 
MODAL MATRIX 

UK : REAL UORK UECTOR OF LENGTH GREATER 
THAN 2*N«N + 3»N 
KPR S PRINT CCNTROL 

=* 0 : PRINT NOTHING 
= l : PRINT EIGENVALUES 
» 2 J PRINT EUERY THING 
FRON : 10-CHARACTER NAME OF CALLING PROGRAM 


i 
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NA.NE : ROM DIMENSIONS 
NAME : 10-CHARACTER NAME FOR MATRIX AA 


OUTPUT ARGUMENTS- 
EIG(N) : 
EE(N.N) : 
EI(N.N) : 


COMPLEX UECTOR OF EIGENUALUE5 OF AA 
MODAL MATRIX MITH ROM DIMENSION NE 
EE INUERSE UITH ROW DIMENSION NE 


CALL GE I GEN ( NWXY , AA . NPDU . 4 . UK . 1 . 2MODAL* . E IG . NWXY . EE . E I . *A~AUCX ) 

RETURN 

END 
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SUBROUTINE PHAG 


>••41 
• • • • « 




THIS ROUTINE CALCULATES AND PRINTS COMPLEX MATRICES 
IN MAGNITUDE AND PHASE FORM. 


CALL STATEMENT ARGUMENTS - 

ARRAY : N X M COMPLEX MATRIX OF ROW DIMENSION NROU 
ITITLE : TUO UCRD LENGTH LITERAL TRRAY CONTAINING THE 
NAME OF ARRAY IN THE MAIN PROGRAM 
AMASE S NRA X 2 ARRAY CONTAINING THE MAGNITUDE AND 

PHASE FORM OF THE ELEMENTS OF THE COLUMNS OF ARRAY 
MATOUT : N X M CGMPLEX MATRIX OF ROW DIMENSION NRA 

CONTAINS ARRAY TRANSFORMED INTO MAGNITUDE AND PHASE 
FORM 


SUBROUTINE FHAG (ARRAY, NROH.N.M. ITITLE. AMASE. MATOUT. NRA) 

COMPLEX ARRAY! NROU. M). MATOUT! NRA. N) 

DIMENSION AMASECNSA.2). ITITLE<2) 

NROU : ROU DIMENSION OF THE INCOMING MATRIX 
NRA : ROU DIMENSION OF THE OUTPUT MATRIX 

PRINT HEADER FOR OUTPUT MATRIX 
BOTH UORDS OF ITITLE CAN BE SET IN A DATA STATEMENT OR 
INDIUIDUALLY USING DIRECT ASSIGNMENT. 

I.E. ITITLE! 1)=10HCONTROLLAB 
ITITLE! 2)=10HILITY====5 

NO ERRORS UILL BE FLAGGED IF ITITLE IS PASSED AS A SINGLE VARIABLE. 
HOUEUER. THE REMAINING UORD UILL EE PRINTED AS JUNK CHARACTERS IN 
THE HEADER. 

URITEIS* 100) 

URITEIG. 110) ITITLE. N.M 
URITEIS, 120) 

PHAG LOOPS THROUGH ARRAY COLUMN BY COLUMN CONUERTING 
EACH ELEMENT INTO MAGNITUDE AND PHASE FORM. 

DO 10 JP-l.N 

REINITIALIZE AMASE 

DO 1 KPM.N 
AMASEOCP. 1)=0.0 
AMASE ! KP » 2 ) =0 . 0 
1 COMTINUE 
DO 5 IP-l.N 

SEPARATE REAL AMD IMAGINARY PARTS OF COMPLEX NUMBERS 
REAL ANO A I MAG ARE ANSI STANOARO FORTRAN FUNCTIONS 
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TERM l =REAL ! ARRAY ( I P . JP ) ) 

TERM2=AIMAG! ARRAY! IP. JP) ) 

AMASEC IP. I )=SORT! TERM1**2+TERM2»*2) 

IF! ABS(TERM1 ) .LT . 1 . OE-15) TERMl=l. OE-15 

ATAN2 RETURMS ARGUMENT IN RADIANS C -PI < THETA < PI ) 

AMASEdP. 2)=ATAN2!TERM2»TERM1 ) 

COMUERSIQM TO DEGREES 

AMASE !IP>2)-57. 29570 » AMASE ( I P . 2 ) 

LOAD IMTO OUTPUT MATRIX 

CMPLX IS ANSI STANDARD REAL TO COMPLEX CONUERSION 

MATCUT ( IP, JP ) =CMPLX ( AMASE ( IP . 1 ) , AMASE ( IP , 2 ) ) 

5 CONTINUE 
10 CONTINUE 

USE COMPLEX MATRIX OUTPUT ROUTINE TO PRINT MATOUT 
HEADER PRINTING SUPRESSEO CLPRINT= -l) 

CALL CMPRNT (NRA. MATOUT, N.M.G, 6. XOUTPUTX.-l) 

100 FORMAT! //v///, 10X.50!*-*)) 

110 FORMAT! 10X.2A10.X MATRIX (X.I2.* BY X.I2.X) X,/, 

■* 15X, ^OUTPUT FORM: ( MAGNITUDE . P iASE(DEG. ) )X) 

120 FORMAT! 10X.50IX-X)) 

RETURN 

END 
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SUBROUTINE RES ID 

****it*»*it»*****«»*«i«*»ti*i«*<**t**t««»*t«t«iiii*tr**»#i»i 

THIS SUBROUTINE DETERMINES THE OUTTER PRODUCT 
OF THE MATRIX OBS AMP A COLUMM OF THE MATRIX COM. 

THIS MATRIX RESI1AT . ^PRESENTS AM IMPULSE RESIDUE 
MATRIX. 

THE ELEMEMTS OF THE RESMAT MATRIX ARE l-IRITTEN TO A TAPE 
SO THAT POST FXECUTIOM ANALYSIS CAN 3E PEPFORMEO. SOME 
EXAMPLES OF THIS MIGHT BE flOPNALIZAT ION OR REORDERING 
OF THE RESIDUES ASSOCIATED WITH PARTICULAR MODES. 


CALL STATEMENT ARGUMENTS - 


COM : NPDU X NC MATRIX OF ROW DIMENSION MUXY 

OBS : NM X fiPDU MATRIX OF RO'.I DIMENSION NZ 

NX : COLUMN DIMENSION OF MATRIX CON 

RESMAT : MATRIX FORMED FROM OUTER PRODUCT OF 

OBS AND A COLUMN OF CON 

I : COLUMN OF CON USED IN DETERMINING RESMAT 


SUBROUTINE RES I D< CON. OBS. NZ» NX. MUXY. RESMAT • I ) 

COMPLEX CON (MUXY. MX). OBS(NZ.MUXY). RESMAT < MUXY . MZ ) 
C 

COMMON ✓DIM-' MU. MC. HD. MM, NP 
C 

NPDU a MP ♦ MO ♦ MU 
C 

DO 10 J=1.MM 
DO £0 K»1»NPDU 

RESMAT(K.J) = OBS(J.K) * COMOC.I) 

30 WRITE! 1. 100) RESMAT ( 1C. J) 

100 FORMAT(aK10.5) 

10 CONTINUE 
RETURN 
END 
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SUBROUTINE SCALE 






IJSJi 


THIS SUBROUTINE SCALES THE NATRICES A. B. C AND D UITH 
A DIACONAL MATRIX T. THIS PROCEDURE IS ANALOCOUS TO 
A CHANCE OF UARIABLE TRANSFORMATION. 


X-DOT = AX ♦ BU ♦ DM 
Y = CX ♦ EU ♦ FU 


APPLY TRANSFORMATION X =* TZ 


2-DOT = (TINU»A*T)Z ♦ (TINU«B)U ♦ CTINU-DH4 
Y * (C»T)X * EU ♦ FU 


SUBROUTINE SCALE (A.B.C.D.NU.NX.NY.NZ.T.TIMU.AS.BS.CS.DS) 
DIMENSION A(NU.NU). B(NU.NX). C(NZ.NU). D(MU.NY) 

DIMENSION AS(NU.NU). BS(NU.NX). CS(NZ.NU). DS(NU.NY) 

DIMENSION TCNU.NU). TINU(MU.NU). URK(SOO) 

COMMON /DIM/ NU. NC. ND. NH» NP 


SUBROUTINE RINU ( TT. TI. UK. NR. N ) 

COMPUTES THE INUERSE OF A REAL MATRIX 
INPUT ARCUMENTS- 

TT(N.N) : REAL MATRIX UITH ROU DIMENSION NR 
UK : REAL UORK UECTOR UITH LENGTH 
GREATER THAN N*N»2»N 
NR : ROU DIMENSION 

OUTPUT ARGUMENTS' 

TKN.N) * INUERSE OF TT UITH ROU DIMENSION NR 


FORM INUERSE OF SCALING MATRIX T 
CALL RIMUIT. TINU.URK.NU.NU) 


SUBROUTINE MP3 ( NA. AR. A. MB. BR. B. NC. CR. CC. C. NP. P ) 

COMPUTE THE MATRIX PRODUCT: P » A«B*C 

A : AR X BR REAL MATRIX UITH ROU DIMENSION NA 

B : BR X CR REAL MATRIX UITH ROU DIMENSION NB 

C I CR X CC REAL MATRIX UITH ROU DIMENSION NC 

P X AR X CC REAL MATRIX UITH ROU DIMENSION NP 
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FORM SCALED DYNAMIC MATRIX AS 
CALL MP3CNW. MU. T.NW. NU. A. MW. NU, NU. TIMU, MW, AS) 


SUBROUT IMF mfq ( AIM. HOUT. MR. NC, MR I. MRO ) 

STORE MATRIX AIM(MR.flC) IN THE MATRIX AOUT(MR.MC) 
WHERE AIM AMD A0U1 HALE ROW 01 HEMS I CMS MR I AMO MRO 


EQUATE SCALED DYNAMIC MATRIX AS TO DYNAMIC 
MATRIX A FOR RE TURN TO MAIM PROGRAM 

CALL MEG(AS.A.MU.MU.NW.MW) 


SUBROUTINE MULT ( A. 0. AB. L. M. M. MA, MB. MAB ) 

COMPUTE MATRIX PRODUCT: AB - A*B 

A : L X M MATRIX WITH ROW DIMENSION NA 

B : M X M MATRIX WITH ROW DIMENSION ND 

A3 : L X M MATRIX WITH ROU DIMENSION NAB 

A. 3. AB REAL MATRICES 


FORM SCALED CONTROL MATRIX BS AND EQUATE IT TO THE 
CONTROL MATRIX B FOR RETURN TO MAIN PROGRAN 

CALL MULTf T. B. BS.NU.NU.NC.MW.NU.NM) 

CALL MEQ(DS.B.ML'.NC.NW.MW) 

FORN SCALED DISTURBANCE MATRIX DS AND EQUATE IT TO 
THE DISTRUBANCE MATRIX D FOR RETURN TO MAIN PROGRAN 

CALL MULT ( T, D. DS.NU.NU. ND.NU.NU.NU) 

CALL M£Q(DS,D.NU.ND.NU.NU) 

FORM SCALED OUTP* * MATRIX CS AND EQUATE IT TO THE 
OUTPUT MATRIX C FUR RETURN TO MAIN PROCRAM 

CALL NULTCC. TINU.CS.NN.NU.NU.N2.NW.MZ) 

CALL MEQiCS.C.NM.NU.NZ.NZ) 

RETURN 

END 
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SUBROUTINE STEP 


THIS SUBROUTINE DETERMINES THE MATRIX UF 
STEP RESIDUES FOR A GIUEN IMPULSE RESIDUE 
MATRIX AND A SPECIFIED SET OF EIGEN- 
UALUES. 

THIS IS DONE BY FORMING A DIAGONAL MATRIX 
OF THE EICENUALUES AND POST - MULTIPLYING ITS 
INUERSE BY THE IMPULSE RESIDUE MATRIX, 


CALL STATEMENT ARGUEMENTS - 

IMP : NPDU X NM MATRIX OF ROW 

DIMENSION NWXY (IMPULSE RESIDUES) 
EIG : NPDU LECTOR OF DIMENSION NWXY 
(UECTOR OF E-UALUES) 

STP : NPDU X NM MATRIX n F ROU 

DIMENSION NWXY (STEP RESIDUES) 
EIGMAT : NPDU X NPDU DIAGONAL MATRIX 
OF ROW DIMENSION NWXY 
(E-UALUES ON DIAGONAL) 

EINU : NPDU X NPDU DIAGONAL MATRIX 
OF ROW DIMENSION NWXY 
(INUERSE OF EIGMAT) 


SUBROUTINE STEP ( IMP, EIG. NWXY, NZ. STP. EIGMAT, EINU) 

CCNPLCX TMP(NUXYfNZ). EIC(NUXY). STP^NUXY, NZ) 
COMPLEX E I GMAT ( NWXY » NWXY ) » EINU (NWXY, NWXY) 

REAL WK(IOQO) 

COMMON /DIM/ NU# NC# ND. NM# NP 

NPDU = NP ♦ ND * NU 

INITIALIZE MATRIX EIGMAT TO ZERO 

DO 5 J=1 • NPDU 
DO G K=l • NPDU 
EIGMAT(J,K) =* 0,0 
CONTINUE 

FORM DIAGONAL MATRIX OF E-UAUJES 

DO 10 I 3 1# NPDU 
EIGMAT (1,1) 3 EIG(I) 


SUBROUTINE CINU ( TT» TI. WK# NR# N ) 
COMPUTE ThE INUERSE OF A COMPLEX MATRIX 


INPUT ARGUMENTS- 
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TT(N*N) : COMPLEX MATRIX UITH ROW DIMENSION NR 
WK : REAL WORK UECTQR WITH LENGTH 
GREATER THAN N»N*2 + 3*N 
NR : ROU DIMENSION 
N : ORDER OF TT AND TI 

OUTPUT ARGUMENTS- 

TICN.N) : COMPLEX INUERSE OF TT WITH 
ROW DIMENSION NR 


INUERT DIAGONAL MATRIX OF E-UALUES 
CALL CINU(EIGMAT»EINU* WK> N WX Y > NPDU ) 


SUBROUTINE MULTCC ( A, B. AB# L# M. N# NA# NB# NAB ) 

COMPUTE MATRIX PRODUCT: AB = A*B 

A : L X M COMPLEX MATRIX UITH ROW DIMENSION NA 

B : M X N COMPLEX MATRIX WITH ROW DIMENSION NB 

AB : L X N COMPLEX MATRIX WITH ROW DIMENSION NAB 


FORM STEP RESIDUE MATRIX 

CALL MULTCC (EINU# IMP, STP» NPDU » NPDU • NM# NWXY * NWXY . NWXY ) 

OUTPUT STEP RESIDUES TO A TAPE 

DO 30 1*1. NM 
DO 20 J= 1# NPDU 
20 WRITE(2» 200) STP(J.I) 

30 CONTINUE 

C 

200 FQRMAT(2F10.5) 

C 

RETURN 

END 
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SUBROUTINE WORK 






THIS SUBROUTINE FIRMS THE COMPLEX MATRICES 
COM. OBS AMD DIST BY MULTIPLYING THE MATICES 
B. C AND D BY A COMPLEX MATRIX M. 

THIS COMPLEX MATRIX MULTIPLICATION REPRESENTS A 
TRANSFORMATION THROUGH THE COMPLEX MATRIX M. WHEN 
THE MATRICES B. C AND D ARE RELATED THROUGH THE 
STATE UARI ABLE EQUATIONS - 

X-DOT - AX + BU + DW 
Y = CX 

SO THAT - 

CON = MINU * B 
OBS = C * M 
DIST = MINU * D 


CALL STATEMENT ARGUMENTS - 


B : NPDU X NC MATRIX OF ROW DIMENSION NWXY 

C : NM X NPDU MATRIX OF ROW DIMENSION NZ 

D : NPDU X ND MATRIX OF ROW DIMENSION NUXY 

M : NPDU X NPDU MARTIX OF ROW DIMENSION NWXY 

MINU : INUERSE OF MATRIX M 

NX : COLUMN DIMENSION OF MATRIX B 

NY : COLUMN DIMENSION OF MATRIX D 

CON : NPDU X NC MATRIX OF ROW DIMENSION NUXY 

OBS : NM X NPDU MATRIX OF ROU DIMENSION NWXY 

DIST : NPDU X ND MATRIX OF ROW DIMENSION NWXY 


SUBROUT I NE WORK (B. C. D.M. MINU.NX.NY. NZ. NWXY . CON . OBS , D I ST ) 
DIMENSION B( NUXY .NX). C(NZ.NUXY). D(NUXY.l) 

COMPLEX M ( NUXY. NWXY ) . M INU ( NWXY . NWXY ) 

COMPLEX CON(NUXY.NX). OBSCNZ. NUXY) . DIST(NUXY, 1) 

COMMON /DIM/ NU. NC. ND. NM. NP 

NPDU = NU + NP + ND 


SUBROTUINE MULTCR ( A. B. AB. L. M. N» NA. NB. NAB ) 

COMPUTE MATRIX PRODUCT? AB = A«B 

A : L X M COMPLEX MATRIX UITH ROW DIMENSION NA 

B ? M X I! REAL MATRIX WITH ROW DIMENSION NB 

AB t L X N COMPLEX MATRIX WITH ROW DIMENSION NAB 
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FORM MATRIX COM 

CALL MULTCRIMINU, B. COM, NPDU, NPDU, NC, NWXY , NUXY, NUXY) 


SUBROUT IME MULTRC ( A, B, AB. L, M, M, NA. MB, MAB j 

COMPUTE MATRIX PRODUCT! AB * A*B 

A ! L X M REAL MATRIX WITH RCll DIMENSION NA 

B ! M X M COMPLEX MATRIX WITH EOU DIMENSION N8 

AB : L X N COMPLEX MATRIX WITH ROW DIMENSION NAB 


FORM MATRIX OBS 

CALL MULTRC (C, M, OBS. NM. N°DU. NPDU, NZ.NUXY.NZ) 

FORM MATIX DIST 

CALL MULTCRIMINU, D, DIST, NPDU, NPDU, 1, NUXY, NUXY .NUXY) 

RETURN 

END 
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